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Abstract 

Neutrino oscillation experiments in the past decades have greatly improved our knowledge on 
neutrinos by measuring the fundamental neutrino parameters. The ongoing and upcoming neutrino 
oscillation experiments are intended to pin down the neutrino mass hierarchy and to discover the 
leptonic CP violation. By means of neutrino oscillograms, we analyze the impact of non-standard 
neutrino interactions on neutrino oscillations in the Earth matter. The standard neutrino oscil- 
lation probabilities may be significantly changed by non-standard interaction parameters, and in 
particular, the CP-violating effects in the energy range = 1 ~ 20 GeV are greatly enhanced. 
In addition, the event rates of muon neutrinos in the proposed huge atmospheric neutrino exper- 
iment, PINGU at the South Pole, have been estimated in the presence of non-standard neutrino 
interactions. It has been found that the PINGU experiment has very good sensitivities to the 
non-standard neutrino interaction parameters. 
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I. INTRODUCTION 



Today, the phenomenon of neutrino oscillations is considered to be the standard and 
leading order mechanism for neutrino flavor transitions, providing strong evidence that neu- 
trinos are massive and lepton flavors are mixed, which leads to physics beyond the Standard 
Model (SM) of particle physics However, although it is widely accepted in the particle 
physics community that this phenomenon stems from a non-trivial structure of leptonic fla- 
vor mixing, the so-called non-standard neutrino interactions (NSIs), which are considered to 
be sub-leading order effects to standard neutrino oscillations, may still affect neutrino flavor 
transitions in a signiflcant way [2j . The concept of NSIs is presently the most studied descrip- 
tion for effects beyond the standard paradigm of neutrino oscillations. In fact, dimension-six 
and higher-order operators exist in various theoretical extensions of the SM, which include 
e.g. seesaw models, i?-parity violating supersymmetric models, left-right symmetric models, 
grand unification theories, and extra dimensions. Basically, all modern extensions could give 
rise to NSIs. Therefore, the investigation of NSIs could reveal additional new physics behind 
neutrino fiavor transitions. In addition, it plays an important complementary role to direct 
searches of physics beyond the SM at colhders such as the LHC. 

For example, NSI effects have previously been studied for the accelerator-based neutrino 
oscillation experiments MINOS and OPERA jsil^], atmospheric neutrino experiments jsj-Q], 
reactor neutrino experiments js iOl, and a future neutrino factory 11-15|. In this work, we 



investigate neutrino flavor transition probabilities based on standard neutrino oscillations 
and NSIs as sub-leading effects. Note that we only consider so-called propagation (or matter) 
NSIs, which are parameterized by different NSI parameters, and not source and detector 
NSIs. Especially, we derive oscillation probabilities for the Vf, — )■ and — )■ channels 
that are important for atmospheric neutrino oscillations, and study the impact of various 
NSI parameters on these probabilities, in both cases of the normal neutrino mass hierarchy 
(NH) and inverted neutrino mass hierarchy (IH). In addition, we estimate the number of 
atmospheric neutrino events in the future PINGU experiment at the South Pole, which has 
been recently shown to have great potential for determining the neutrino mass hierarchy 



16Hl9|. Very good sensitivities to the NSI parameters are expected at this experiment. 

The remaining part of our paper is organized as follows. In Sec. II, we will review the 
formalism of three-flavor neutrino oscillations in the presence of standard and non-standard 
matter effects, and present the parameter mappings between the leptonic mixing matrix in 
vacuum and that in matter. The analytical approximate formulas of oscillation probabilities 
are derived for the appearance channel Vf. — )■ z/^ and disappearance channel Vy, — )■ z/^. Sec. Ill 
is devoted to numerical analyses of the NSI effects on the neutrino oscillation probabilities 
with and without intrinsic CP violation. Using neutrino oscillograms, we further point 
out the important regions in the plane of zenith angles and neutrino energies, where the 
NSI effects are most signiflcant. Taking into account the NSIs, we estimate the event rate 
of atmospheric muon neutrinos at the PINGU detector in Sec. IV. It turns out that the 
PINGU has very good sensitivities to the NSI parameters. Finally, we conclude in Sec. V. 
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II. NEUTRINO OSCILLATIONS WITH NSIS 



In general, NSIs can be present for neutrino production, propagation, and detection. We 
will concentrate on the non-standard matter effects for neutrino propagation, which should 
be relevant for atmospheric and long-baseline neutrino oscillation experiments in particular. 
In this section, we recall the general formulation of three-flavor neutrino oscillations in the 
presence of standard and non-standard matter effects. In this case, neutrino flavor transitions 
are determined by the effective Hamiltonian 



(1) 



where the vacuum Hamiltonian is given by 
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the standard matter potential is 

Uraix) = Vcc 

with Vcc = V^GpNeix), and the non-standard matter potential is parametrized by 
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where are real (for a = f3) or complex (for a ^ (3) constants, i.e., the so-called NSI 
parameters. Note that Gp is the Fermi constant and N^[x) is the electron number density 
in matter. 

The exact oscillation probability is given by PaS = |5'/3a(a;, Xo)^, where Si3a{x,XQ) is the 
evolution matrix satisfying the Schrodinger-like equation. 



i-\u{x))^n{xMx)) , 

and can be obtained by solving the above equation as 

S{x,Xo) = exp — i / 1-L{x') dx' 



(5) 



(6) 



Since the standard matter potential "Hm is invariant under any rotation in the 2-3 plane, 
it is sometimes convenient to work in a new flavor basis {i^e-i^2i^iV — ^23(^^6' '^m' '^t)^ 
with C/23 = 023-^5- Note that the leptonic mixing matrix can be parametrized as U — 



3 



^23-^(5^13-^1^12! where O^^ denotes the rotation in the i-j plane with a rotation angle 
and = diag(l, 1, e'^) with 6 being the leptonic Dirac CP-violating phase. If neutrinos are 
Majorana particles, there will be two additional CP-violating phases in U. However, these 
two leptonic Majorana CP-violating phases are irrelevant for neutrino oscillations both in 
vacuum and in matter. 

Now, in this basis, the standard Hamiltonian can be written as 



/o 

Am|i 
^0 








^12^13 + 
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where V-sDix) = UI^'Hsb{x)U23 with 7/30(2;) = "Ho + 'Hra{x), while the non-standard matter 
potential turns out to be 



Is 
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with the modified NSI parameters 
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Hence, in this basis, the effective Hamiltonian is 
/ 



Uix] = A 
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where A31 = Aml^/{2E), a = Aml^/Aml;^ and A = Vcc/A3i. Additionally, note that 

Cj • have been defined. The evolution matrix S{x, Xq) in this 



s,„- = sin6',-- and cos6',,,- 

'-J 



basis is related to S{x, Xq) in the flavor basis via the unitary transformation S{x, Xq) = 
1/2^3 {x,Xo)U23. The oscillation probabilities of antineutrinos can be obtained through the 
replacements A — )■ —A and 6 — > —6. 



A. Parameter Mappings 



20 



and 



21 



we use a perturbation method to derive the effective 



Now, following Refs. 

neutrino masses and leptonic mixing matrix [7™ in matter. Note that current neutrino 



oscillation data indicate a 



0.03. Therefore, we keep the terms of the same order 



4 



a and 5^3, and ignore all other higher-order contributions, such as as^3, a^, and asl^,- Thus, 
we write V. = M ■ A^^ and introduce M = A^(°) + A^^^^: 



13'-13 
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where A^*^*^) corresponds exactly to the standard Hamiltonian "Hgj-, in the two-flavor limit with 
a = 0, while A^'-^^ incorporates the corrections from a and the NSI parameters. Obviously, 
AA^^^ can be diagonalized by a rotation in the 1-3 plane, i.e.. 



/ cos sin^^gX 



V- 



sin6']^3 cos 6^^ J 



(12) 



with the efi^ective mixing angle 613 given by 



tan 26* 



sin 26 



13 



13 



COS 26*^^3 — A 



(13) 



In the following, we choose 9^^ to be deflned in the flrst quadrant ^. Namel y, G [0,7r/4] 
for A < cos26'i3 and 6^^ G [tt/A,tt/2] for A > cos26-^^^. Therefore, we obtain 

C - cos 2^13 + A 
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sin' 9 
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C + COS 26*13 - ^ 
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(14) 



with C = a/ (cos 29y^ — Ay + sin^ '^^u- Two other useful relations can readily be derived 
from Eq. (14) 

cos 2^13 = , (15) 



sin 2^ 



sin 2^ 



13 



13 
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implying sin 6*13 = sin6'i3/(l — A) at the leading order of sin 6*13. In the limit of a = and 
in the absence of NSIs, ^13 is just the effective mixing angle in matter. 
Furthermore, the eigenvalues to zeroth order are given by 

1 
2 

. 



(0) 



1 + A-C 



(0) 
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A + C 
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^ Note that it is also possible to define 613 to be in [0,7r/4] by properly arranging the eigenvalues and the 
corresponding eigenvectors, as shown in Ref. [22 1. 
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and the eigenvalues to first order are A 
forward to sliow 



(1) 



(1) 



a COS^ ^135^2 + ^ 



Mf^ witli = t/(o)^A^(i)f/(o). It is straigfit- 
cos^ ^i3£^ee ~ sin 2^^3Re(eg^) + sin^ ^136^ 



.(1) 



.(1) 
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a sin^ ^13^12 + A 



sin^ 9i3e„„ + sin 29...Re(e„^) + cos^ Oise. 



The corrections to tlie eigenvectors are given by 
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wfiere t/^*'^^ and 11^^^ stand for the column vectors of the matrices U^^^ and f/*^°^ , respectively. 
The effective neutrino masses in matter are determined by rhf = ml + Aml^lXf'^ + A-^^], 
while the effective leptonic mixing matrix is U"" = U^^ + U^^^ . In the absence of NSIs, 
there are two resonances, i.e., A = a and A = cos2^^3. For neutrino energies E > 1 GeV 
and matter densities p = 3 g/cm^ in the Earth crust, we have A > a, so only the resonance 
A = cos 26^^ is relevant 22| . Now, after some lengthy computations, we find that 



= sin + {sin 2^13 [asj, + A{e,, - 5, J] + 2 A [cos 2^i3Re(5,J + ilm(5,j] } 
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q;C23 sin26'^2 sin6']^3 
l + ^ + C 



(19) 



which reproduce the well-known results in the limit of vanishing NSI parameters 
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Our results in Eq. (19) differ from those in Ref. 23|] by including higher-order corrections 



23|. 



from a and Sy^. The parameter mapping is not valid in the resonance region, where the 
perturbation theory breaks down. Furthermore, the mixing matrix elements in matter could 
be divergent, so we will calculate the oscillation probabilities that should be well-defined in 

the mixing 



22 



general, and particularly in the resonance region. As pointed out in Ref. 
angle Qyi can be arbitrary in the limit of a = and in the absence of the NSIs, so the mapping 
for f/^ in Eq. (19) cannot be taken seriously. For neutrino energies E > \ GeV, onlythe 
effective mixing angle ^13, which is determined by [/^ in the standard parametrization [l|, is 
crucially important for neutrino oscillations in the Earth matter. Note that the true effective 
mixing angle ^13 differs from ^13 in the contributions from a and the NSI parameters. 



B. Oscillation Probabilities 



In practice, the oscillation probabilities can be computed using perturbation theory based 
on small quantities, e.g., the smallest mixing angle 6'i3, the ratio of the two mass-squared 
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differences a, and the NSI parameters Eafs- Following Ref. [2l|, one can explicitly decompose 
the effective Hamiltonian as 1-L{x) = TIqIx) + 7ii{x), where 'Hq(x) is the zeroth order in the 
small parameters and 'Hi{x) includes the higher-order contributions. Then, to first order, 
the evolution matrix is approximately given by 

S{x,xo) c:^ So{x,xo) - iSo{x,xo) / [So{x',xo)~'^Hi{x')So{x',xo)]dx' , (20) 

Jxo 

where the zeroth-order evolution matrix Sq{x,xq) is determined by T-Lq{x). For constant 



matter density, the oscillation probability for the Ue — )■ i^fj, channel is found to be [sl. Il2| 



^ea^ — ^eu ~ 4S13S23C23 (|£e/i|c23Cx — l^er | ■S23Caj) 



sm sin h sin 

2 2 2 



+ 8S13S23 [\ee^,\S23Cx + |£er|C23Ca;J ^ _ ^ Siu ^ 

/, , , , , AA A (l-A)A . , 

+ 8S13S23C23 \\^ep\c2-iSx - \eeT\s2-iS^^) sm — sm — sm ^ , (21) 

where A = A-^^^L denotes the oscillation phase in vacuum, and §13 = Si3/(1 — A) is just sin 6^13 
in the leading order. In addition, we have defined e^^ = I^Q^je''^"'' (for a/3 = e/x, er), x = 
+ 6, and u = + 6. Note that we have neglected the terms proportional to a e^^, and 
Pg^^ stands for the transition probability without the NSIs, i.e, P^J;^ ~ ^§^3^23 sin^ ^^~2^^ ■ 
The approximate formulas of neutrino oscillation probabilities to the second order of a and 



A . 2(1-^)A 



Si3 can be found in Ref. [2l|. At leading order, only the NSI parameters and appear in 
the transition probability Pf^^ [sl, Esl- Hence, we will concentrate on these two parameters in 
following numerical analysis. Furthermore, the CP-violating terms in the last line of Eq. (21) 
related to NSI parameters are not suppressed by the ratio of two mass-squared differences 
a, compared to the standard case. Thus, even when the standard CP violation is not visible 
in an experimental setup, one may expect observable CP-violating effects stemming from 
the catalysis of NSIs. 

Next, for the z/^ — )• i/^ channel, we have 

m 

A 



P^^'-Pl^ - k^rlc^,. (sL23^AsinA + 4s2x23cL23^sin2 



2 



+ (k/,;,! - |£..) 1 5^^23^2x23 sin A- 2v4sin2^^ , (22) 

where e^r = |e^r|e"'^''^, ^2x23 = sin 2^23 (^2x23 = cos2^^23 have been defined. Note that 



only the NSI parameters e^^, e^,-, and appear in the survival probability Pj^^^ [sj. Since 
the current experimental bound on is very stringent, the dominant NSI effects should 



come from e^,- and Ej 



III. NUMERICAL ANALYSIS 

In order to illustrate the NSI effects on neutrino propagation in the Earth, we calculate 
numerically the effective mixing angle 613 and the oscillation probabilities. In our numerical 



7 




I 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



-0.3 -0.2 -0.1 0.1 0.2 0.3 




1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



FIG. 1: Dependence of the effective mixing angle sin^ ^13 on the NSI parameters and neutrino 
energy, where a constant matter density profile p = 4.5 g/cm^ has been assumed (i.e., the Earth 



mantle density) and the dotted line cor resp onds to sin^ ^13 
values of neutrino parameters from Ref. 2g] have been used. 



sin 01 3. In addition, the best-fit 



computations, we assume the PREM model of Earth matter density 2J], and express the 
baseline in terms of the zenith angle as L = —2Rcos6z, where R = 6371 km is the Earth 
radius and 6^ = tt — h with h being the nadir angle. Note that —1 < cos^^ < —0.84 
corresponds to the trajectories crossing both the mantle and core of the Earth, while —0.84 < 
cos ^2 < to those crossing only the Earth mantle. On the other hand, there already exist 
restrictive experimental constraints on the NSI parameters in realistic models 25|. However, 



in Ref. 26|, the model-independent upper bounds on the matter NSI parameters have been 



found to be much larger than the model-dependent ones: 



/ keel < 4.2 


ke^l < 0.33 


kerl 


< 3.0 \ 




\e^^\ < 0.068 


k/ir 


< 0.33 


\ 




l^TT 1 


< 21 J 



(23) 



Therefore, in the following discussions, we just ignore e^^, which receives the most stringent 
constraint. For the other matter NSI parameters, we will take a conservative value {eapl =0.1 
for illustration ^. 



A. Effective Mixing Angle in Matter 

First of all, it may be interesting to show how the standard and non-standard matter 
effects modify the effective neutrino mixing angles 6ij in matter. As mentioned in the 
previous section, we will focus on 613, which is relevant for neutrino energies E > 1 GeV. To 



^ Recently, the MINOS experiment has constrained the NSI parameter to the range —0.20 < e^,- < 0.07 at 



the 90% confidence level 27[, in the framework of two-flavor neutrino oscillations. 
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examine the dependence of ^13 on the NSI parameters, we consider two specific examples, 
where only one relevant NSI parameter is switched on in each case and all the CP-violating 
phases (i.e., 6 and (pa/s) are set to zero: 

• Est 7^ 0. With the help of Eqs. (9) and (19), we derive 

■?7{ (7 - cos 26*13 + A cos 26*13 -A, . /^.s 

Sin^ ^13 = + ^ A Sm 2^i3C23£er , (24) 

26 6'^ 

where the small terms asi3 and e^^ have been neglected. It is worthwhile to mention 
that if the resonance condition A = cos 26*13 is satisfied, we have ^13 = 7r/4, which is 
independent of the NSI parameter in the leading order approximation. This can 
be well understood in the framework of two-fiavor neutrino oscillations in matter with 
NSIs 0, 29|, where one observes that the off-diagonal term in Hnsi cannot modify the 



resonance condition. The result in the case of Eefj. 7^ can be obtained by replacing 
C23£eT with 8236 in Eq (24), and the difference between these two cases can be at- 
tributed to a non- maximal ^23- In Fig- El we have calculated sin^ 6*13 using the exact 
formulas of parameter mappings |23|. From the left panel, we can clearly observe that 
the resonance condition is essentially unchanged by Ser, and ^13 = ^13 is achieved when 
the standard matter effects are cancelled by the NSI effects. However, the resonance 
is in fact shifted by higher-order corrections. 

£fiT 7^ 0. In a similar way, we obtain 

. 2~ C- COS 2^13 + A sin^ 2^13 , . _„ 

sm ^13 = ^ ^^^Asm2^23£^r , (25) 

where the small terms asi^ and have been omitted. It is now evident that the 
standard resonance condition A = cos 26*13 does not lead to sin^6*i3 = 1/2. Namely, 
the resonance has been shifted by the NSI parameter, which has also been pointed 
out in Ref. j29| in the framework of two-neutrino oscillations. For a fixed value of 
^13, if Efj^r becomes larger, the neutrino energy has to increase in order to balance the 
negative contribution from e^r- In the right panel of Fig. [H we have shown sin^ 6*13 
by using the exact mapping formulas. The main features can be well described by the 
approximate formula in Eq. (25). Note that, compared to Eer, the e^r correction to 
^13 is milder since it is further suppressed by sin 26*13. Similarly, one can also consider 
the impact of See, £^^1, and Err on the effective mixing angles. 



B. Neutrino Oscillograms: Standard Case 



The matter effects on neutrino propagation in the Earth can be perfectly illustrated 
through the so-called neutrino oscillograms. In order to compare between standard and 
non-standard matter effects, we first briefiy summarize the general features of the standard 
neutrino oscillograms, which have been systematically studied in Refs. |30|, |31 
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FIG. 2: Standard neutrino oscillograms without NSIs (i.e., £ai3 = 0) in the appearance channel: 
= P{i'e — >• f/i) for neutrino oscillations in the case of normal neutrino mass hierarchy (left 
panel) and P^^ = P{i'e — ^ for antineutrino oscillations in the case of inverted neutrino mass 
hierarchy (right panel). 
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FIG. 3: Standard neutrino oscillograms without NSIs (i.e., Eap = 0) in the disappearance channel: 
P{v^ — )• v^) for neutrino oscillations in the case of normal neutrino mass hierarchy (left 
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panel) and in the case of inverted neutrino mass hierarchy (right panel). 



In Figs. [2] and |3l we have reproduced the neutrino oscillograms in the Ve — ^ and 
^// — ^ t'/i channels, respectively. In our numerical calculations, the latest global-fit data on 

have 



28 



leptonic mixing angles and the neutrino mass-squared differences given in Ref. 
been used. Since the Dirac CP- violating phase has not been experimentally constrained, we 
simply take 5 = in the calculations. Some general comments on the neutrino oscillograms 



are in order 30 



In the left panel of Fig. |51 the oscillation probability = P^Ue ^ ^n) is shown in 
the NH case. The resonance in the mantle appears around cos^^ ~ —0.75 and E-^ ^ 
6 GeV, where the resonance energy is essentially determined hj A = cos 26i3 while 
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the corresponding baseline Lr or the zenith angle 6^ by the requirement of maximal 
oscillation phase A = tt. Note that we can safely ignore the effects of Amgi for neutrino 
energies E > 1 GeV. The ridges in the core region are caused by the parametric 
resonances, receiving both contributions from the mantle and core oscillation phases 



30 



. Since the resonance takes place in the neutrino sector, the oscillation probability 
of z/g — )■ Ufj, is suppressed by matter effects and the oscillogram is almost empty. In 
the right panel of Fig. [2|, we have calculated the antineutrino oscillation probability 
= P{'^e ~^ '^fi) in the IH case. Now that the resonances occur in this case and 
dominate the contributions to oscillation probabilities, the similarity between P^^ for 
NH and P^^ for IH is evident. The oscillation probability Pf^ in the NH case is highly 
suppressed, as P^^ in the IH case. 

• In Fig. [31 we have given the survival probabilities = P(z^^ — )■ i^fj,) for both NH (left 
panel) and IH (right panel). The probability P^^ receives the dominant contribution 
from the vacuum oscillation due to 6*23 and Am|i, and is significantly affected by the 
1-3 mixing through 6*13 only in the resonance regions. In fact, the latter effect reduces 
P^^, as indicated in the left panel of Fig. [31 In the absence of resonance, as in the 
IH case, P^j^ is basically given by the vacuum oscillation probability as shown in the 
right panel. The antineutrino survival probabilities P^^ = Pi'^n ^n) in NH and IH 
cases are similar to P^^ in IH and NH cases, respectively. 

The detailed study of the maxima and minima in neutrino oscillation probabilities and the 



general conditions for resonances in the two-flavor approximation can be found in Ref. 



30 



while those for the three-flavor oscillations in Ref. 3l| . In our discussions, since the neutrino 



energies are always above 1 GeV, the three-flavor corrections should be negligible. 



C. Neutrino Oscillograms: NSI Effects 

Now we proceed to discuss the NSI effects on the standard neutrino oscillations in the 
Earth by using neutrino oscillograms. To quantify the NSI effects, we consider the difference 
between the standard (i.e., Sap = 0) and non-standard (i.e., eai3 7^ 0) neutrino oscillograms, 
namely AP^^ = P^^^ - Pf^ and AP^^ = Pj^^^ - P^°. In Fig. [H we numerically calculate 
APe^ and ^P^i^i in the NH case, and only one NSI parameter is switched on in each plot. 

First, we consider the difference APe^ in the appearance channel, for which the approx- 
imate formula in the case of constant matter density can be obtained from Eq. (21) and 
the numerical results are shown in the left column of Fig. [H Two comments are in order: 
(1) The contributions from NSI parameters to APg^ are proportional to S13 in the limit 
y4 — )■ 1, i.e., around the resonance. In the deep core, i.e., cosOz < —0.84, the parametric 
resonances dominate, so the NSI effects are not necessarily suppressed. (2) The difference 
APe^ is independent of e^r in leading order. For this reason, as shown in the lower plot, 
APe^ is vanishing everywhere except in the core, if only e^j^r 7^ is assumed. One common 
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FIG. 4: Differences between tlie standard and non-standard neutrino oscillograms in the z/e — )• z/^ 
channel APe/. = Pe^^ - (left column) and in the i/^ channel AP^^ = P^^^i _ pSD ^-^igj^^ 

column), where the normal neutrino mass hierarchy and (5 = are assumed, and ee/^ = 0.1 (upper 
row), Eer = 0.1 (middle row), and e^r = 0.1 (lower row) are taken for illustration. 



feature of all three plots is that significant NSI effects are lying in the core region, which is 
only accessible in the atmospheric neutrino experiments. 

Second, we turn to the difference AP^^ in the disappearance channel, for which the 
approximate formula in the case of constant matter density can be obtained from Eq. (22) 
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and the numerical results are shown in the right column of Fig. HI As implied by Eq. (22), 
the effects induced by Eefj. and Eer can only arise from higher-order corrections of S13 and a, 
so they are insignificant as in the upper and middle plots. The most interesting observation 
is that I AP^^I ~ 1 could be achieved, particularly in the core region. Switching off both 
and Err in Eq- (22), we have 

AP^^ ^-\E^r\ sin 2^23 (^AA sin^ 2^23 sin A + 4 cos^ 2^23 sin^ , (26) 

where = is assumed. Furthermore, note that the second term in the parentheses on 
the right-hand side of Eq. (26) is always positive, while the first term can be either positive 
or negative, depending on the oscillation phase A. If A = {2k + l)vr/2 with k being an 
integer, the contribution from the second term is negligible, since cos^ 26'23 ~ 0.03 for the 
best-fit value of ^23 = 40°. Therefore, AP^^ oc sin A should show an oscillatory behavior, 
as given in the lower plot. It is worthwhile to note that A = {2k + l)7r/2 corresponds to 
the oscillation minima of the standard survival probability P^^ ■ For example, comparing 
the left plot of Fig. [3] with the lower plot in the right column of Fig. SI one can observe 
the huge difference in the region along the diagonal line. If A = /cvr holds, the second term 
dominates over the first one, and AP^^ follow the same oscillatory structure as the leading 
vacuum oscillation term in P,?,?. 



D. Non-Standard CP Violation 

Finally, we show that the NSIs may lead to the enhancement of CP-violating effects in 



neutrino oscillations |32j. In Fig. [5l the probability differences APf,^{5) = P^f^^{S) —P^f^^{0) 
induced by the leptonic CP violating phase 6 are calculated for 6 = 7i/2. To signify the 
CP-violating effects due to 6, we set all the NSI phase parameters to zero, i.e., (pap = 0. As 
mentioned before, the CP-violating effects in the standard case come from the interference 
between two different oscillation frequencies, and thus these effects are relatively small in 
both cases of NH and IH. However, the probability differences are enhanced greatly when 
the NSI parameters are switched on, in particular in the NH case. 

The left and right plots in the first row correspond to APe^{6) = P^^{6) — P^^{0) in 
the NH and IH cases, respectively, since the NSI parameters are taken to be zero. One can 



observe that the impact of S is insignificant. See Ref. [33[ for more discussions about the 
CP-violating effects in the standard three-neutrino oscillations in matter. Even if the NSI 
parameters are nonzero, the CP-violating effects are negligible in the IH shown in 

the middle and lower plots in the right column. This is because the transition probability 
P^if^ itself is suppressed rather than enhanced by matter effects. However, as in the NH 
case, |APe/i(^)| for S = 7i/2 could be 12% in the mantle, while as large as 30% in the core, as 
illustrated in the middle and lower plots in the left column. At the same time, without the 
NSI effects, it is vanishingly small in both mantle and core regions. Therefore, the future 
long-baseline oscillation experiments, where neutrino beams traverse the Earth mantle, will 
have excellent sensitivities to NSI-enhanced CP-violating effects. 
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FIG. 5: Oscillograms for the probability differences APe^{6) = Pijr (S) - Pil^^O) with 6 = tt/2 in 
the normal neutrino mass hierarchy case (left column) and in the inverted neutrino mass hierarchy 
case (right column). The two plots in the first row correspond to APe^('^) with vanishing NSI 
parameters, i.e., APe/i('^) 



P!^iS)-P!^iO) with<5 = 7r/2. 



The remarkable difference between standard and non-standard CP-violating effects be- 
comes clear, if we look at the last line of Eq. (21). More explicitly, APen{6) is proportional 
to asi3 in the standard case, while to l^e^ilsis or |£:er|si3 in the non-standard case. For 
£ai3 = 0.1 and a ~ 0.03, the standard CP violation is suppressed by one order of magnitude. 
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IV. EVENT RATES AT PINGU 



In previous sections, we have explored the general features of NSI effects on neutrino 
oscillations in the Earth matter at the level of oscillation probabilities, which should have 
important implications for atmospheric neutrino experiments. In order to study the NSI 
effects in a realistic experiment, we have to calculate the number of neutrino events in 
small bins of energies AE and zenith angles A cos 6^ for a terrestrial detector. To this end, 
we take a multi-megaton scale ice Cherenkov detector for example, such as the proposed 
PINGU detector at the South Pole 



34 



The main idea of this proposal is to make the 
Deep Core of the IceCube detector denser to lower the energy threshold down to a few GeV, 
implying possible precision measurements of atmospheric neutrinos. Such an experimental 
setup would be rather interesting in view of its great physics potential for the determination 
of neutrino mass hierarchy, and the oscillation parameters 16l4l9|. We refer the readers to 
Ref. 3J] for a detailed description of the PINGU experiment. 

Atmospheric neutrinos interact with the nucleons via charged-current interactions in the 
detector and produce energetic charged leptons, which are radiating Cherenkov photons 
when propagating in ice. The Cherenkov photons will be captured by the dense strings of 
photomultipliers. The charged muons /i^ leave clear and long tracks in the detector, so 
the experimental resolution to the direction of muon neutrinos is much better than that 
of electron neutrinos. Now, we make a rough estimate of the distribution of i/^-like events 
collected by the PINGU detector for one year running. Explicitly, the number of z/^-like 
events in the i-th zenith-angle bin and j-th energy bin is given by 



d cos 6, 



dE 



+ 



d$, 



pNSI I 

nil ~U 



Ej 

d^. 



d^„ 



pNSI I 
;/// ~r 



d cos OzdE 



NSI 



P 

d cos ME 



pNSI 

d COS ME 



a 



cc 

PN 



pVesiE) 



(27) 



d cos ME '''' 

where <l>y^ {^u^,) denotes the neutrino z/q, (antineutrino z/q,) fluxes, A^a is the Avogadro's 
number, P^^^ (E'a^^) stands for the neutrino (antineutrino) oscillation probabilities^ with 
NSI effects, and the effective volume of PINGU with 20 strings is parametrized as 

E y'-' 
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pVes{E) = 14.6 Mt X 



log 



10 



GeV J 



(28) 



Furthermore, we adopt the following simple parametrization of the deep inelastic u-N and 
p-N scattering cross sections 
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a^S{E) = 7.30 X 10"^^ cm^ 



a 



cc 

PN 



(E) 



3.77 X 10" 



cm^ 



VGeV 

(— 

VGeV 



(29) 



The atmospheric electron and muon neutrino fluxes are taken from Ref. [35j], where the 
neutrino fluxes have been calculated using a hybrid method for 1 GeV < E < 10^ GeV and 
tabulated in small bins of neutrino energy and zenith angle. 
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008 6*^ COS 02 

FIG. 6: The number of f^-like events in the E — cos9z plane, collected by the PINGU detector for 
one year, where the normal neutrino mass hierarchy is assumed. The upper left plot corresponds 
to the case of standard neutrino oscillations, i.e., all the NSI parameters are vanishing. 



The event distribution in the case of standard neutrino oscillations is shown in the upper 
left panel of Fig. [6l where the bin widths AE = 1 GeV and Acos6'^ = 0.05 are used 
and NH is assumed. Note that the highest number of events can reach ~ 800, which is 



in agreement with the estimation in Ref. [16[. It is evident that the event distribution is 
mainly determined by the z/^ — > survival probability, which is given in the left panel of 
Fig. [3l In the presence of NSIs, the event distribution is distorted, as shown in the other 
plots of Fig. [6l In particular, in the case of e^r 7^ 0, the number of events at high energy 
bins is increased remarkably. This is because the z/^ — z/^ channel in Eq. (!27|) dominates the 
contributions to the observed z/^-like events, and the NSI parameter e^,- significantly modifies 
the survival probability in this region. The latter has already been pointed out in Sec. Ill 
C and illustrated in Fig. |3l Therefore, the PINGU detector has a better sensitivity to 
than Ee/i and Eer- Note also that the total number of events is approximately unchanged. 

To quantify the significance of NSI effects, we define an SD-NSI asymmetry as the dif- 
ference between the number of z/^-like events in the standard case (i.e., = 0) and that 
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A for specific values of NSI parameters. In the cases of Ee^ and Eer 7^ 0, the large 
asymmetry appears in the low energy regions E ~ [5, 10] GeV, and is separated by the lines 
of zero asymmetry. For E^fj, = 0.1 or Eer = 0.1, the maximal asymmetry can be as large as 



7, which is comparable to the estimated sensitivity to the neutrino mass hierarchy [16 
As expected, in the case of a non- vanishing Ef^r, the SD-NSI asymmetry is much more 
stunning (e.g., A ~ 50), indicating a great discovery potential on this NSI parameter. We 
have also made a comparison between NH and IH in Fig. [TJ where one can observe that the 
asymmetry in some bins is reduced in the case of ^e^i 7^ and Eer 7^ for IH, while there 
is no significant difference in the case of e:^,- 7^ 0. Therefore, the experimental sensitivity 
to the NSI parameters, in particular E^r, at the PINGU detector, is almost independent of 
neutrino mass hierarchy. 

Note that the real sensitivity to NSI effects will actually be much lower because of sys- 
tematics, the smearing effects in the reconstruction of neutrino energies and zenith angles, 
the uncertainties in the matter density profile and the other neutrino oscillation parameters 



16l |. In addition, the Ur interactions lead to an important background to the events. 
Nevertheless, a detailed simulation on the probe of NSIs at PINGU is meaningful and will 
be elaborated elsewhere. 



V. SUMMARY 

Now we are entering a new era of precision measurements of neutrino parameters, includ- 
ing three leptonic mixing angles and two neutrino mass-squared differences. The ongoing 
and forthcoming neutrino experiments are expected to pin down neutrino mass hierarchy 
and to discover CP violation in the lepton sector. As the precisions of oscillation experiments 
gradually improve and more data are accumulated, it is promising to discover or constrain 
the new physics effects beyond the standard paradigm of neutrino oscillations. 

One of the most widely studied scenarios is the non-standard neutrino interactions. In 
this paper, we have considered the NSI effects on neutrino oscillations in the Earth matter 
by using neutrino oscillograms. First, we derive the mapping formulas between the lepton 
mixing matrix in vacuum and that in matter beyond the leading order approximation. In 
particular, the NSI effects on the effective mixing angle 61^ have been discussed in some 
detail. Then, the NSI effects on the neutrino oscillograms of the Earth are investigated. The 
most significant difference between the standard and non-standard oscillograms appears in 
the z/^ — )■ z/^ survival probability and in the case of 7^ 0. In addition, the CP-violating 
effects in neutrino oscillations can be enhanced by the NSI effects, even if the CP-violating 
phases in the NSI parameters are switched off. Finally, the NSI effects in the PINGU 
experiment are explored. We calculate the event rate of atmospheric muon neutrinos at 
PINGU, and demonstrate that the future huge atmospheric neutrino experiments should 
have very good sensitivities to the NSI parameters, in particular E^r- However, a more 
sophisticated simulation of the NSI effects at PINGU, including the systematics and other 
uncertainties, is needed to make a final conclusion. 
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In addition to atmospheric neutrino experiments, such as PINGU, the ongoing and up- 
coming long-basehne neutrino osciUation experiments, which are intended for the determi- 
nation of neutrino mass hierarchy and leptonic CP violation, are also sensitive to the NSI 
effects. In general, the future neutrino oscillation data will soon lead us either to the discov- 
ery of new effects beyond the standard oscillation scenario, or to more restrictive constraints 
on new physics parameters. 
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